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ABSTRACT 

We present new ATCA observations at 3.3 mm of 27 young stellar objects in the p— Oph young cluster. 25 of these 
sources have been detected. We analyze the sub-millimeter and millimeter SED for a subsample of 17 isolated class II 
protoplanetary disks and derive constraints on the grain growth and total dust mass in the disk outer regions. All the 
disks in our sample show a mm slope of the SED which is significantly shallower than the one observed for the ISM at 
these long wavelengths. This indicates that 1) class II disks in Ophiuchus host grains grown to mm/cm-sizes in their 
outer regions, 2) formation of mm/cm-sized pebbles is a fast process and 3) a mechanism halting or slowing down the 
inward radial drift of solid particles is required to explain the data. These findings are consistent with previous results 
in other star forming regions. We compare the dust properties of this sample with those of a uniformly selected sample 
in Taurus- Auriga and find no statistical evidence of any difference in terms of grain growth between the two regions. 
Finally, in our sample the mm slope of the SED is not found to correlate with indicators of grain growth to micron 
sizes in the surface layers of the inner disk. 
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1. Introduction 

Circumstellar disks around pre-main sequence (PMS) stars 
are the common birth places of planetary systems. In or- 
der to build up large bodies such as giant planets a huge 
growth of solid particles by more than 12 orders of mag- 
nitude in size has to occur starting from submicron-sized 
microscopic dust grains as those typically found in the in- 
terstellar medium (ISM; Mathis et al. H977|) . The first stages 
of this process of grain growth are characterized by the dy- 
namical interaction between gas and dust, leading to col- 
lisions between the solid particles and finally coagulation 
(see Beckwith et al. |2"0"0UI Dominik et al. |2T)0"71 Natta et 
al.[MZD. 

Evidence for the presence of micron-sized dust grains 
in protoplanetary disks has been provided by the inspec- 
tion of ISO and Spitzer spectra for the silicate feature at 
about 10 fim (e.g. Bouwman et al. 200TJ Van Boeckel et 
al. 120031 Kessler-Silacci et al. [2006[) . However infrared ob- 
servations can only probe the uppermost surface layers of 
the disk, while to investigate the properties of dust in the 
disk midplane, where planet formation is expected to take 
place, observations at longer wavelengths are needed. If the 
long-wave emission is optically thin the spectral index of 
the SED at these wavelengths can be related to the spec- 
tral index of the dust opacity j3 (ftdust ^ an d P- values 
lower than ~ 1 are naturally interpreted in terms of grain 
growth (DraineHDnH. Beckwith & Sargent (fTWT) were the 
first to show that T Tauri Stars (TTS) have SEDs at sub- 
millimeter wavelengths which are typically shallower than 
what found for the ISM, suggesting that dust grains as large 
as at least 1 mm are present in the disk midplane. However 
these single-dish observations did not have a good enough 
angular resolution to spatially resolve the disks, and the 



same data could in principle be explained also by optically 
thick emission from compact disks with unprocessed, ISM- 
like grains. 

Wilner et al. (pOOHf and Testi et al. (|2"005|l resolved 
the disks around TW Hya and CQ Tau PMS stars 
at 7 mm with the Very Large Array (VLA), confirm- 
ing that the long-wave emission from these disks reveals 
the presence of cm-sized pebbles in the disk midplane . 
More recently Rodmann et al. (Rodmann et al. (2006)1 
and Lommen et al. (|2007[) found evidence of grain growth 
to ~ mm/cm-sized pebbles for about ten T Tauri stars 
in the Taurus-Auriga, Chamaeleon and Lupus star form- 
ing regions (SFRs). Lommen et al. (201fJ| have combined 
new and literature data for disks from five different SFRs 
(Taurus- Auriga, Lupus, Chamaeleon, Corona Australis and 
Serpens) and reported a tentative correlation between the 
mm slope of the SED and the strength of the 10-/im sili- 
cate feature, possibly suggesting that inner and outer disk 
evolve simultaneously in terms of dust grain growth. 

So far, most of the observations carried out at long- 
wavelengths have targeted the brightest sources. Ricci et 
al. (|2010l hereafter R10) have started to probe fainter disks 
with new sensitive data at 3 mm obtained with the Plateau 
de Bure Interferometer in Taurus and analyzed a sample of 
21 isolated class II Young Stellar Objects (YSOs) in this 
SFR. For all the disks in this sample they found neither 
evidence for time evolution of dust grain growth nor sig- 
nificant relations between dust properties and stellar ones 
were found. 

In this paper we perform the same analysis presented 
in R10, but on a sample of disks in another SFR, i.e. 
p— Ophiuchi. In Sect. 2 we present new 3 mm data for 
27 p— Oph YSOs obtained with the Australia Telescope 
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Compact Array (ATCA!]]) and its new Compact Array 
Broadband Backend (CABB). In Sect. 3 we describe the 
properties of the sub-sample used for the analysis. The 
method adopted for deriving the dust properties, namely 
grain growth and dust mass, is the same used in RIO and 
it is described in Sect. 4. The results of the analysis are 
described in Sect. 5, whereas Sect. 6 summarizes the main 
findings of our work. 

2. New ATCA observations 

We observed 27 YSOs in the p— Oph star forming region 
at 3.2944 mm with ATCA and the new CABB digital filter 
bank. Targets were chosen by being class II (most of them) 
or "flat spectrum" YSOs, with a detected flux at 1.3mm 
higher than 20 mJy in order to have good chances for de- 
tection at ~ 3.3 mm with a few hours at most of integration 
time on-source. 

The observations were carried out between the end of 
April and the beginning of July 2009, with the ATCA ar- 
ray in the H168 and H75 configurations respectively. The 
obtained FWHM of the synthesized beam of - 3 - 7" did 
not allow us to spatially resolve any of our sources, and so 
a 2D-gaussian model was used to estimate the flux density 
in the continuum at 3.3 mm. To obtain the best sensitiv- 
ity in the continuum we set the correlator to cover the full 
2-GHz CABB effective bandwidth on all the five antennas 
available at ~ 3 mm. Each source has been observed one or 
more times, with a total integration time on-source ranging 
between about 40 minutes and 3 hours. 

The complex gain calibration was done on the calibra- 
tors 1622-297 and 1730-130, and science targets were typ- 
ically observed for 10 minutes and then spaced out with 
2 minutes on the gain calibrator. The absolute calibration 
was done with Uranus except for a few times in which 1921- 
293 was used and then its flux was always cross-calibrated 
with Uranus the day after. In general, the baselines that 
resolved out Uranus were not used for flux calibration^. In 
this work, we assume a la— uncertainty on the calibrated 
flux of about 25%. 

The results of the observations are reported in Table [T] 
In particular we have detected 25 out of 27 targeted YSOs. 
The detected sources have fluxes ranging from about 1.5 
and 48.8 mJy at 3.3 mm. For the two undetected sources, 
IRS 33 and GY 284, we obtained an upper limit for the 
flux at 3.3 mm of 1.2 mJy. These sources have been previ- 
ously detected by Motte et al. (fTi?9"8"|> using the IRAM 30m 
single-dish telescope at 1.3 mm. Although their reported 
fluxes are relatively large (105 ± 20 mJy and 130 ± 10 mJy 
for IRS 33 and GY 284 respectively) the emission for both 
the sources appears to be spatially resolved by a beam with 
a FWHM of about 11 arcsec. This indicates that a signif- 
icant fraction of the collected flux at 1.3 mm comes from 
an extended envelope surrounding the circumstellar disk. 
In order to get a reliable estimate for the flux at ~ 1 mm 
from the circumstellar disk only, interferometric observa- 
tions at ~ 1 mm are needed to filter out the contaminating 
emission from the envelope. Since these are not available in 

1 The Australia Telescope Compact Array is part of the 
Australia Telescope which is funded by the Commonwealth 
of Australia for operation as a National Facility managed by 
CSIRO. 

2 For the data reduction we used the MIRIAD package. 



the literature we do not include IRS 33 and GY 284 in our 
sample discussed in the next section. 

3. Sample 

In this section we describe some properties of the sample 
considered for our analysifl In Sect. 13.11 we outline the 
selection criteria adopted to select our sample, Sect. 13.21 
describes the sample completeness level, and finally in Sect. 
13.31 we derive the stellar parameters. 

3.1. Selection criteria 

We used the same selection criteria adopted in R10: 

1. class II YSOs as catalogued in Andrews & 
Williams (|2007a[) from the infrared SED to avoid 
contamination of (sub-)mm fluxes by a residual 
envelope; 

2. central stars need to be well characterized through 
optical-NIR spectroscopic/photometric data to obtain 
self-consistent disk SED models; 

3. YSOs with at least one detection in the 0.45 < A < 1.3 
mm spectral rangeQ other than being observed through 
the new ATCA observations described in Sect. 2, to have 
a good sampling of the (sub-)mm SED; 

4. no evidence of stellar companions with projected phys- 
ical separation between 5 and 500 AU to avoid tidal 
interactions that may alter the outer disk structure; 
adopting for all the YSOs in our sample the current 
estimate of ~ 130 pc for the distance to the p— Oph 
star forming region (Wilking et al. 120081 Lombardi 
et al. I2008|) . the range in projected phyical separa- 
tion translates into a range in angular separation of 
- 0.05 - 4". 

Among the 25 YSOs detected with ATCA at 3.3 mm, 
17 satisfy these selection criteria and they constitute the 
sample of our analysis. They are listed in Tabled 

3.2. Completeness 

The two histograms in Fig. [1] show the distribution 
of the "isolated", i.e. satisfying the selection criterion 
(4) in Se ct 1511 class II YSOs from the Andrews & 
Williams (2007a) catalogue. The histogram on the left side 
shows that our sample comprises all the isolated class II 
YSOs with Fi.3 mm > 75 mJy, while for the fainter objects 
the completeness level reduces to 22% (8/36). 

In terms of the stellar properties, our sample includes 
53% (16/30) of the isolated PMS stars with stellar types 
equal to or earlier than Ml, correspondent to a stellar mass 
larger than ^0.5 M for a ~ 1 Myr old PMS star according 

3 Note that the final sample used in the analysis does not 
comprise all the targets of the new ATCA observations described 
in Sect.d 

4 Note that in R10 a slightly narrower spectral range of 0.45 
< A < 0.85 mm was chosen. We adopted a broader spectral 
range to include in our sample the sources YLW 16c and WSB 
52, for which no observations at A < 1.3 mm have been carried 
out so far. The different uncertainties for the spectral slopes of 
these sources, due to a shorter spectral leverage, are properly 
taken into account in the analysis. 
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to the Palla & Stahler (fTOM)) models^ (see Fig.0). Our sam- 
ple is much less complete for PMS stars with later spectral 
types (or smaller stellar mass): only 7% of the known PMS 
stars with a mass smaller than ~ 0.4 M Q is in our sample 
(1/15), and this has to be considered as an upper limit since 
the Andrews & Williams (|2007a|l catalogue does not com- 
prise class II YSOs with spectral types later than M5. This 
is due to the sensitivity limits of the current (sub-) mm facil- 
ities to detect emission from circumstellar material around 
very-low mass PMS stars. 



3.3. Stellar properties 

To constrain the stellar properties we used the same method 
as described in RIO. We estimated the stellar effective tem- 
peratures by converting the adopted spectral types (from 
Andrews & Williams I2007al and references therein^) with 
the dwarf temperature scale of Schmidt-Kalcr (1982) for 
spectral types earlier than MO and the intermediate tem- 
perature scale of Luhman (1999) for spectral types equal to 
or later than MO. 

Stellar luminosities were computed from the 2MASS J- 
band flux (Cutri et al. 12003) after calculating the extinction 
of each object by dereddening the J — H and H~K S colors 
to the locus observed for Classical T Tauri stars (Meyer et 
al. I1997|) . and adopting the Cardelli et al. (|1989p extinction 
law with Ry — 4.2, which is appropriate for p— OplQ. 

Luminosities and effective temperatures were con- 
verted into stellar masses and ages by using the Palla & 
Stahler (1999) models of PMS stars, as done in RIO (see 
Fig. [2]) . According to these models the ranges spanned by 
our sample go from about 0.2 M Q (WSB 60) to 1.9 M Q 
(SR 21) in mass and from about 0.7 Myr (GSS 26, EL 24, 
YLW 16c) to 6.6 Myr (SR 21) in age. 

The stellar parameters are reported in Table [5J 
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Fig. 2. H-R diagram for the sources of our sample. The 
dashed and solid lines represent the isochrones and evolu- 
tionary tracks respectively from the Palla & Stahler (1999) 
PMS evolutionary models. In this diagram the evolutionary 
tracks start from an age of 0.1 Myr. The uncertainties, not 
shown in the plot, are typically ~ 0.08 — 0.13 in logL* and 
~ 100 — 150 K in temperature. 



4.1. Disk structure 

In order to characterize a model of the disk estimates for 
some stellar physical quantities (bolometric luminosity L*, 
effective temperature T e ff and mass M+), plus some infor- 
mation on the disk structure and on the the dust opacity 
are needed. As for the stellar parameters we used the val- 
ues listed in Table [21 Regarding the disk structure the only 
relevant parameters for our analysis are the disk outer ra- 
dius i? ut and the parameters S<j U st,i and p which define a 
power-law surface density for the dust component: 



4. Analysis 

In order to constrain the dust properties in the disk outer 
regions we fitted the sub-mm/mm SED with two-layer (i.e. 
disk surface and midplane) models of flared disks heated by 
the radiation of the central star (Chiang & Goldreich 1997, 
Dullemond et al. I2001|) . In Sect. [O] and 03] we describe 
the parameters which are needed to define a two-layer disk 
model. In Sect. 4.3 we then discuss which are the physical 
quantities of these models that can be constrained by our 
analysis. 



5 Note that this completeness level is very close to the one 
obtained in R10 in the Taurus-Auriga star forming region, i.e. 
~ 58% of the isolated class II YSOs with stellar masses larger 
than ~ 0.4 M©. 

6 For only two cases in our sample, i.e. EL 20 and IRS 41, the 
spectral types are not available in Andrews & Williams ( 2007a). 
For these PMS stars we adopted the stellar types from Natta 
et al. (2006J which are based on near infrared broad band pho- 
tometry. 

7 For IRS 41, for which the 2MASS J-band flux is not 
available, we adopted the luminosities as derived by Natta et 
al. (|2006l) . after correcting them by a multiplicative factor (130 
pc/150 pc) 2 to account for the different adopted distance of 150 
pc in the Natta et al. paper. 



Sdust(-R) - S dust ,i f - Au j , (1) 

radially truncated at i? ut- For the disk outer radii we 
adopted the intervals listed in Column (3) of Table [31 For 
the seven disks in our sample which have been mapped by 
high angular resolution observations at sub-mm/mm wave- 
lengths presented in Andrews & Williams (12007b!) we con- 
sidered intervals roughly centered on their best-fit values. 
Three sources (EL 20, WSB 52, DoAr 44) have been ob- 
served and spatially resolvedby recent observations carried 
out by Andrews et al. (|2009|l with the SMA at 0.87 mm: 
these disks show evidence of an outer radius larger than 
100 AU, but an accurate estimate of the outer radius i? ut 
could not be derived, since in order to fit the interfero- 
metric visibilities they adopted the self-similar solution for 
a viscous disk (Lynden-Bell & Pringle I1974p rather than 
a truncated power-law disk. For the seven sources in our 
sample which have not been mapped yet we adopted a 
fiducial interval for _R out of 100 — 300 AU, which com- 
prises nearly all the disks outer radii as found for the 24 
Taurus and Ophiucus circumstellar disks in Andrews & 
Williams (2QQ7b). Based on the results of this same work 
we adopt for all the sources in our sample an interval for 
the surface density power-law index p which goes from 0.5 
to 1.5. 
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Fig. 1. Histograms higlighting some properties of our selected sample. In both the histograms our sample is represented by 
black columns, while the total columns (black+white) include all the class II YSOs from the Andrews & Williams (|2007ap 
catalogue with no evidence of stellar companions in the 0.05 — 4.0" interval in angular separation, from which our sample 
has been selected (see Sect. I3.1| ). Left) Distribution of the fluxes at 1.3 mm, including upper limits; Right) Distribution 
of stellar spectral type (see Sect. 



An important thing to keep in mind here is that with 
these values for the outer radius of the disk the dust emis- 
sion at the long wavelengths considered in this paper turns 
out to be dominated by the optically thin outer disk re- 
gions. This has two important consequences for our discus- 
sion. The first one is that constraints on dust properties 
such as dust grain sizes and dust mass can be derived from 
the continuum emission at sub-mm/mm wavelengths. If a 
disk is much more compact than the ones which have been 
commonly mapped so far, i.e. if i? ut 20 — 30 AU, then its 
emission would be dominated by the denser inner regions 
which are optically thick even at these long wavelengths, 
and no information on the dust properties could be inferred 
by the observed continuum. The second consequence is that 
the disk inclination is not a relevant parameter for our anal- 
ysis, except only for the case of a nearly edge-on disk that 
however would make the central PMS star invisible in the 
optical. 

4.2. Dust opacity 

To calculate the dust opacity we adopted the same dust 
grain model taken in R10, i.e. porous composite spheri- 
cal grains made of astronomical silicates (optical constants 
from Weingartner & Draine l^OOip . carbonaceous materials 
(Zubko et al. I1996P and water ices (Warren [1984| with frac- 
tional abundances from a simplification of the model used 
in Pollack et al. (1994) and a volume fraction for vacuum 
of w 30%. In the disk surface and midplane we consider 
a dust grain population with a grain size number density 
throughout all the disk 

n(a) cx a~ q (2) 

between a m ; n and a max (a here is the grain radius). In this 
paper we adopt two different grain size distributions in the 
disk surface and midplane, since both observations and the- 
ory predict larger grains in the midplane. In particular, al- 
though we consider in both the surface and midplane a 
small value of a m i n ~ 0.1 /jm, for the maximum grain size 
cimax in the disk surface we consider a value only slightly 



larger than a m iJ3i whereas a max has been left free to vary 
in the disk midplane. As done in R10 we consider in this 
work four possible values for the q— parameter: 2.5, 3.0, 3.5, 
4.0. 

4.3. Method 

After having set the disk outer radius i? ut and the power- 
law index p of the surface density profile, the two-layer mod- 
els can be used to fit the sub-mm/mm SED of circumstellar 
disks to constrain dust properties in the disk midplane^., in 
particular the spectral index of the dust opacity /3 between 
two mm- wavelengths (1 and 3 mm in the case of this paper) 
and the product Md us t x Kimm between the mass in dust and 
the dust opacity at 1 mm. Except for very low values of /3, 
which cannot be explained by large values of the power-law 
index q of the grain size number-density (see discussion in 
Sect. 5.1), every /3- value can be reproduced by different (g, 
ttmax) couples. Furthermore, at a fixed q, the precise value 
of the maximum grain size a max correspondent to a certain 
derived (3 depends strongly on the model that one adopts 
for the dust. The only robust conclusion which is valid for 
all the reasonable models of dust analyzed so far is that /3- 
values lower than the value found for the ISM (/3tsm ~ 1.7) 
can be obtained only with dust populations in which grains 
as large as at least ~ 1 mm are present (see e.g. Natta et 
al. I2007p . For these reasons in the rest of the paper we will 
rarely refer to the maximum grain size a max , whereas we 
will more frequently use /3 as our proxy for grain growth. 

In order to derive an estimate for Ki mm and thus con- 
strain the dust mass Mdust, one has to adopt a certain 
model for the dust grain which provide a family of func- 
tions Ki mm (/3) labeled with the g-parameter. The Ki mm (/3) 



We also consider the same value of q for both the surface and 
the midplane. However, contrary to the dust in the midplane 
where q is a relevant parameter, in the disk surface the adopted 
q— value is practically non influential since a max is only slightly 
larger than a m in 

9 Note that this is true only in the case in which the dust 
emission at long wavelengths comes mostly from the optically 
thin disk outer regions, as described at the end of Sect. 14. ll 
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Table 1. Summary of the 3.3mm ATCA observations. 



Object 


a 

(J2000) 


8 

(J2000) 


F v 
(mJy) 


rms 
(mJy) 


Comments 15 


Ref. 


Other names 


AS 205A 


16:11:31.4 


-18-38-26 


27.2 


1.2 


triple (1.5", sb) 


1 


HBC 254 


SR 4 


16-25-56 1 


-24:20:48.3 


4.4 


0.4 






AS 206 


GSS 26 


16-26-10 3 


-24:20:54.9 


24.2 


0.9 








EL 20 


16-26-18 9 


-24:28:20.2 


7.3 


0.5 






VSSG 1 


LFAM 1 


16:26:21.7 


-24:22:50.8 


17.5 


0.6 


flat spectrum 


2 




DoAr 24E 


16:26:23.4 


-24:21:00.7 


8.3 


0.4 


binary (2.0") 


1 


GSS 31 


DoAr 25 


16:26:23.7 


-24:43:14.1 


25.0 


0.6 






WSB 29 


EL 24 


16:26:24.1 


-24:16:14.0 


48.8 


0.7 






WSB 31 


EL 27 


16:26:45.0 


-24:23:08.2 


38.7 


0.6 






GSS 39 


WL 18 


16:26:49.0 


-24:38:25.7 


3.1 


0.5 


binary (3.6") 


3 


GY 129 


SR 24S 


16-26-58 5 


-24:45:37.1 


26.6 


0.8 


flat spectrum, triple (5.2", 0.2") 


2,4,5 


HBC 262 


SR 21 


16:27:10.2 


-24:19:12.9 


4.2 


0.4 






EL 30 


IRS 33 a 


16:27:14.5 


-24:26:46.1 


< 1.2 


0.4 






GY 236 


IRS 41 


16-27-19 3 


-24:28:44.4 


6.2 


1.2 






WL 3 


CRBR 85 


16:27:24.7 


-24:41:03.2 


1.5 


0.3 


envelope 


6 




YLW 16c 


16:27:26.5 


-24:39:23.4 


6.5 


0.4 






GY 262 


GY 284 a 


16:27:30.8 


-24:24:56.0 


< 1.2 


0.4 








IRS 49 


16-27-38 3 


-24-36-58 8 


4.4 


0.4 






GY 308 


DoAr 33 


16:27:39.0 


-23:58:19.1 


3.7 


0.3 






WSB 53 


WSB 52 


16-27-39 5 


-24-39-15 9 


10.2 


0.5 






GY 314 


IRS 51 


16:27:39.8 


-24:43:15.0 


12.7 


0.8 


binary (1.6") 


3 


GY 315 


WSB 60 


16:28:16.5 


-24:36:58.0 


15.3 


0.5 






YLW 58 


SR 13 


16:28:45.3 


-24:28:19.2 


10.0 


0.6 


binary (0.4") 


3 


HBC 266 


DoAr 44 


16:31:33.5 


-24:27:37.7 


10.4 


0.5 






HBC 268 


RNO 90 


16:34:09.2 


-15:48:16.9 


7.6 


0.5 






HBC 649 


Wa Oph 6 


16:48:45.6 


-14:16:36.0 


10.3 


0.4 






HBC 653 


AS 209 


16:49:15.3 


-14:22:08.6 


17.5 


0.5 






HBC 270 



a) For this undetected source the (a, S) coordinates are from 2MASS (Cutri et al. I2003[> ; 

b) Reason why the source has not been considered in the analysis (see Sect. 13. I'll : enclosed in parenthesis there are the angular- 
separation^) between the stellar companions ("sb"= spectroscopic companion); 

c) References: (1) McCa be et al. (f2006)l : An drews fc Williams (20073); (3) Ratzka et al. ([2005jl : (4) Reipurth & Zinnecker lfT993)l : 
(5) Simon et al. (1995); (6) McClure et al. pOlO) . 



functions considered in this paper have been obtained us- 
ing the dust model presented in Sect. 14.21 and they are the 
same shown in Figure 3 in R10. 

Figure [3] reports the best fit flared disk model over- 
plotted to the sub-mm/mm data for each disk in our sam- 
plc0- Since as explained in Sect. 4.1 for each disk we have 
adopted, instead of a single value, an interval of possible 
values for i? ut and p, the uncertainty on these parameters 
translates into an uncertainty on the quantities derived by 

10 Note that for a few disks (SR 4, EL 20, RNO 90) the flux 
at ~ 1 mm falls below the model line. Although this could be 
due to problems in the observations, another possible reason is 
that the 3 mm-flux for these sources is contaminated by free-free 
emission. If this was the case, the derived /3- values for these disks 
would be only lower-limits. Observations at longer wavelengths, 
where free-free dominates the emission, are needed to constrain 
its possible contribution at 3 mm. 



fitting the sub-mm/mm SED, i.e. /3 and Md us t X K lmm . 
Adding this contribution to the uncertainties in the obser- 
vational data, the total absolute uncertainties are approxi- 
mately 0.4 for (3 and a factor of « 3 — 4 for Mdust x i^imm- 

5. Results 

The results of the SED-fitting procedure are listed in 
Table 3. Here we discuss them in terms of dust grain growth 
(Sect. 5.1) and mass in dust (Sect. 5.2) for our sample of 
protoplanetary disks. 

5.1. Grain growth 

Information on the level of dust grain growth in the outer 
regions of protoplanetary disks comes from the analysis of 
the spectral index of the dust opacity at (sub-) millimeter 
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Table 2. Stellar properties of the considered sample. 



Object 


ST 


T c ff 




R, 


M* 


Age 






(K) 


(Le) 


(Re) 


(M e ) 


(Myr) 


SR 4 


K5 


4350 


2.17 


2.58 


1.14 


1.07 


GSS 26 


MO 


3850 


1.39 


2.64 


0.56 


0.54 


EL 20 


MO 


3850 


0.93 


2.16 


0.62 


1.07 


DoAr 25 


K5 


4350 


1.43 


2.10 


1.12 


2.09 


EL 24 


K6 


4205 


2.58 


3.01 


0.96 


0.58 


EL 27 


MO 


3850 


0.78 


1.98 


0.58 


1.23 


SR 21 


G3 


5830 


11.38 


3.29 


1.97 


2.21 


IRS 41 


K7 


4060 


1.61 


2.55 


0.80 


0.76 


YLW 16c 


Ml 


3705 


1.11 


2.55 


0.48 


0.54 


IRS 49 


MO 


3850 


1.02 


2.26 


0.64 


1.00 


DoAr 33 


K4 


4590 


1.81 


2.12 


1.44 


2.88 


WSB 52 


K5 


4350 


0.95 


1.71 


1.04 


4.17 


WSB 60 


M4 


3270 


0.23 


1.49 


0.20 


0.93 


DoAr 44 


K3 


4730 


1.55 


1.85 


1.29 


5.13 


RNO 90 


G5 


5770 


10.24 


3.19 


1.87 


2.32 


Wa Oph 6 


K6 


4205 


2.32 


2.86 


0.98 


0.71 


AS 209 


K5 


4350 


2.11 


2.55 


1.18 


1.23 



wavelengths /3, which reflects the spectral index of the disk 
SED at these long wavelengths a. In particular for a com- 
pletely optically thin disk in the Rayleigh- Jeans regime 
f3 = a — 2, whereas if emission from the optically thick 
disk inner regions (i.e. R < 20 — 30 AU) and deviations of 
the emitted spectrum from the Rayleigh- Jeans regime are 
taken into account, as done for our analysis, ft > a — 2. In 
the limit case of a completely optically thick disk even at 
these long wavelengths the SED spectral index a becomes 
independent on /? and so no information on grain growth 
would be obtainable. In Columns (5) and (6) of Table 3 
the constrained values of a and /3 between 1 and 3 mm are 
reportecO 

Figure 4 shows the SED spectral index ai_3 mm plotted 
against the observed flux at 3.3 mm for all the sources in 
our sample. The spanned range in ai_3 mm is 1.8 — 2.9, and 
there is no clear evidence of any correlation between the two 
plotted quantities. The fact that all the disks in our sample 
show a (sub-)mm spectral index which is shallower than 
the one found for the ISM, i.e. a < aisM, brings evidence 
of grain growth from an initial ISM-like dust population for 
all the disks in our sampleP^I. 

In Figure 5 we report the histogram of the derived dis- 
tribution of /3-values. All the disks show (3 < /3ism ~ 1<7, 



11 Note that the difference a — ft turns out to be in the range 
1.6—1.9. This (small) discrepancy from the value of 2 is primarly 
due to the low temperature of the outer disk midplane and the 
consequent deviation from the Rayleigh- Jeans regime of the mm- 
wave emission. 

12 This sentence is strictly valid for the ten disks in our sample 
which have been mapped so far; for the seven disks which have 
not been mapped yet this sentence is valid only if the underlying 
assumption that their spatial extension is not so small that most 
of the mm emission is optically thick (7? ut J; 20 — 30 AU) holds 
true. 
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Fig. 4. Spectral index a between 1 and 3 mm from the 
model best fits shown in Figure 3 plotted against the ob- 
served flux at 3.3 mm. Open blue points reprensent the 
spatially resolved disks, whereas the filled red points are 
for the unmapped ones. The dashed vertical line shows the 
mean 3tr-sensitivity of our ATCA observations, i.e. about 
1.8 mJy. 



indicative of dust grain growth to at least mm-sizes. Note 
that for 9 out of the 17 disks j3 < 0.8 (considering a ler- 
uncertainty on j3 of about 0.4) and for them the data are 
not consistent with the MRN-value of 3.5 for the power- 
law index q of the dust grain size distribution in the ISM 
(Mathis et al. 11977ft : for these disks q— values as low as 
2.5 — 3 are needed (see discussion in R10). 
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Fig. 3. Sub-mm/mm SED fits for the 17 sources in our sample with the two-layer flared disk models (solid lines). The 
errorbars in the plots take into account an uncertainty of 10% on the absolute flux scale at wavelengths shorter than 3 
mm, and of 25% at 3 mm. The fitting values of the spectral index a and of the dust opacity spectral index /3 between 1 
and 3 mm are indicated in the bottom left corner of each plot. The value of the adopted surface density power-law index 
p is 1, the adopted value for the outer disk radius is the central value of the i? ou t — interval listed in Table |31 whereas the 
value for the power-law index of the grain size number-density q is 3 except for four disks (GSS 26, YLW 16c, IRS 49, 
WSB 52) for which models with q = 3 cannot explain the data (see Table |3]). For these sources models with a q— value 
of 2.5 arc shown. 



The mean value is < j3 >= 0.46 ± 0.08, which is consis- discussed in RIO (< /3 >Tau= 0.6±0.06). In order to make a 
tent within la with the one obtained for the Taurus sample statistical comparison of the dust grain growth between the 
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Fig. 5. Distribution of the dust opacity spectral index /3 
for the disks in our sample. In blue the values for the spa- 
tially resolved disks are indicated, whereas the red is for 
the unmapped ones. The average /3 value for all the sources 
in our sample and the value of ~ 1.7 for the ISM dust are 
indicated as dashed vertical lines. 

Taurus sample and the one presented in this paper in the 
p— Oph star forming region we performed a two-sample KS 
test. The probability that the two samples have f3— values 
drawn from the same distribution is « 22%. The hypoth- 
esis that the samples in Taurus and Ophiuchus have the 
same /3— distribution cannot thus be rejected at the 95% of 
confidence level. 

In Figure 6 we plot the (sub-)mm spectral index of the 
dust opacity j3 against the stellar age for YSOs in different 
evolutionary stages: the 38 class II disks in Ophiuchus and 
Taurus presented here and in RIO respectively, and a sam- 
ple of 15 less evolved class YSO^f] in Taurus, Perseus, and 
isolated from the PROSAC survey ( Joergensen et al. I2007P 
and from Kwon et al. (2009). Contrary to the class II disks, 
for all the class objects the values of /3 have been obtained 
using the approximated ft = a — 2 relation (see caption of 
Figure 6). As described before, this relation gives only a 
lower-limit for /3, and this probably explains why for many 
of these sources the derived (3 is negative. A more sophisti- 
cated analysis is needed to get more robust estimates of (3 
by taking into account self-consistently deviations from the 
above relation as due to departures from the Rayleigh- Jeans 
regime of the emission (expecially for these cold sources) 
and to marginally optically thick emission typically associ- 
ated to the compact structure forming the disk. However, 
the low values of f3 obtained for nearly all these sources 
appear to show evidence for dust grain growth to ~ mm- 
sizes (see Figure 3 in RIO) already in the earliest stages 
of star formation. Ormel et al. (j2009|) have recently inves- 
tigated the effects of dust coagulation and fragmentation 
onto the dust size distribution in molecular cloud cores. 
They found that grain sizes close to ~ 1 mm can be formed 
if cloud lifetimes are not restricted to free-fall times but 
rather support mechanisms like e.g. ambipolar diffusion are 
present and if frcczc-out of ice has commenced. According 
to their simulations ice-coated grains can grow to sizes of 
~ 0.3 — 8 mm in one ambipolar diffusion timescale at den- 

13 Note that for these objects there are no robust age-estimates. 



sities of n = 10 5 — 10 7 cm -3 , which are typical of the inner 
regions of molecular cores. Dust grain growth to ~ mm-sizes 
can thus be a process accompanying the very first phases 
of star formation. 




0.1 1 10 
toga (Myr) 

Fig. 6. Dust opacity spectral index j3 versus age: relation- 
ship between the dust opacity spectral index /3 and the 
estimated stellar age obtained as in Sect. 3.3 for the class 
II YSOs. Note that for the class YSOs (green circles and 
black stars) the position in the x-axis of the plot is only 
indicative, since no robust age estimates can be derived for 
these sources. Blue squares represent the sample of class II 
YSOs presented in this paper, red asterisks are the Taurus- 
Auriga class II disks from R10, green circles are the class 
YSOs from the PROSAC survey (Joergensen et al. lMJTj) . 
black stars are class YSOs from Kwon et al. (12009I) . The 
(3- value for the ISM dust is indicated as a dashed horizontal 
line. For all the class YSOs the values of j3 have been ob- 
tained using the f3 = a — 2 relation, where a is the spectral 
index between the CARMA total fluxes at 1.3 and 2.7 mm 
for the Kwon et al. sample, and between the fluxes collected 
at 0.85 and 1.3 mm with the SMA array at baselines longer 
than 40 kA for the PROSAC sample. This criterion was 
chosen by the authors to minimize the contribution from 
the extended envelope to the total emission. 

Figure 6 shows also that there is no relation between the 
dust opacity spectral index and the stellar age for the class 
II disks: grains as large as ~ 1 mm appears to be present in 
the outer regions of disks throughout all the class II evolu- 
tionary stage. This is in contrast with the short timescales 
of inward radial drift expected for ~ mm/cm-sized grains 
in the outer disk as a consequence of the dust interaction 
with the gas component. In order to explain the retention 
of large dust grains in these outer regions some mechanisms 
which may halt the drift of solid particles, e.g. local pressure 
maxima due to turbulent vortices or spiral density waves, 
have been invoked. Birnstiel et al. (|2010b|) have compared 
the observed fluxes at millimeter wavelengths for the disks 
samples described in R10 and in this paper with predic- 
tions of dust evolution models accounting for coagulation 
and fragmentation. They showed that, if radial drift of solid 
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particles is completely suppressed, a grain size distribution 
at the steady-state (due to a balance between coagulation 
and fragmentation) can explain the mm-wave emission of 
the brightest disks. The observed flux of the fainter disks 
are instead typically overpredicted even by more than one 
order of magnitude. These discrepancies may be explained 
by considering in the disk models a dust reduction due to 
radial drift at a reduced rateP^I or during an earlier evo- 
lutionary time or due to efficient conversion of dust into 
larger, unseen bodies (see Birnstiel et al. I2010bl for more 
details). Observations of these faint disks can thus help us 
to determine which mechanisms play a major role for the 
dynamics and evolution of large grains in the outer regions 
of disks. 
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Fig. 7. SED spectral index between 1 and 3 mm versus in- 
dicators of grain growth from the 10-/im silicate feature. 
Left) SED slope between 1 and 3 mm vs the ratio between 
the 10- /^m line flux and the continuum as derived from 
Spitzer/IRS observations by Furlan et al. (2006) for the 
sample of Taurus- Auriga class II YSOs described in R10. 
Right) SED slope between 1 and 3 mm vs the ratio between 
the flux at 11.3 and 9.8 ^m as derived from Spitzer/IRS ob- 
servations by McClure et al. (2010); the plotted points are 
all the class II YSOs described in this paper with a value 
of F11.3/F9.8 reported in McClure et al.. In the two plots 
the arrows in the lower left corners show toward which di- 
rection the values of the plotted quantities are indicative of 
the presence of largest grains. 

Figure 7 shows the spectral index between 1 and 3 mm 
in Taurus (left) and Ophiuchus (right) plotted against two 
indicators of dust processing from the silicate feature ob- 
served at about 10 fim (see e.g. Kessler-Silacci et al. 2006). 
Data for the 10 /jm silicate feature for the two samples in 
Taurus and Ophiuchus come from the literature and re- 
fer to different indicators (see caption of Figure 7) since 
we could not find the same indicator for both the regions. 
However Lommen et al. (|2010l) showed that these two indi- 
cators correlate well and so they can be both used to probe 
the growth of grains from interstellar, submicron sizes to 



The radial drift, other than decreasing the amount of dust in 
the disk and thus decreasing the flux at millimeter wavelengths, 
is more efficient for mm/cm-sized pebbles than for smaller grains 
in the outer disk. For this reason, if one wants to explain the low 
values of the mm-spectral indeces only a reduced rate of radial 
drift (from that expected theoretically) can be invoked. 



sizes of several microns in the disk surface layers. In par- 
ticular low values of the ratio between the 10 /iin line flux 
and the continuum, (F10 — -FcontV-Pcont, and large values of 
the ratio between the fluxes at 11.3 and 9.8 /im, Fn. 3/Fg.s, 
are interpreted in terms of grain growth to micron sizes 
(Bouwman et al. 120011 Kessler-Silacci et al. 2006). Very re- 
cently Lommen et al. (2010) have found a tentative correla- 
tion between these quantities for a sample of about 30 TTS 
and Herbig-Ae/Be systems spread over five different star 
forming regions, including Taurus-Auriga. Since the mm 
slope of the SED probes grains in the disk outer regions 
whereas the 10-/im silicate feature is sensitive to grains in 
the inner regions, the tentative correlation may indicate a 
parallel evolution of the inner and outer disk in terms of 
dust grain growth. However, in this work we do not find 
any correlation neither for the sample in Taurus nor for the 
one in Ophiuchus. Note that in the case of Ophiuchus the 
sample with literature data for the 10- /im silicate feature is 
limited to eight disks only. In Taurus our sample comprises 
nine of the eleven disks considered by Lommen et al., and 
for two of them (AA Tau and GM Aur) our derived values 
of ai_3 mm are not consistent with and lower than the values 
used by Lommen et al0- Note also that the absence of such 
a correlation would not be too surprising since the physical 
mechanisms which are responsible for the observed values 
of the grain growth indicators are different: the presence of 
mm/cm-sized pebbles in the midplane of the outer disk is 
mainly due to coagulation processes and mechanisms which 
trap these pebbles in the outer disk, whereas the presence 
of /im-sized grains in the surface layers of the inner disk is 
probably regulated by fragmentation of larger solid parti- 
cles and a balance between settling and turbulence mixing 
which keeps these relatively small grains in the uppermost 
layers of the disk. Further observations with more sensitive 
telescopes in the future are needed to provide necessary the 
necessary statistics to better investigate the possible rela- 
tion between dust grains in the outer and inner disk regions 
and to possibly constrain the processes of radial mixing and 
vertical settling. 

5.2. Disk mass 

Table 3 lists the derived dust masses for our sample of disks 
detected at 3 mm. As already discussed in R10, the in- 
ferred dust mass depends, at a given chemical composition 
and porosity for the dust grain, on the assumed value of 
the power-law index q of the grain size distribution. This 
dependence, which is due to the different millimeter dust 
opacities obtained for different values of q (see R10), is par- 
ticularly strong for relatively low values of (3. For example, 
in the case of WSB 60, with /? rj 0.3, is larger than 

^dust by a factor of about 10. 

Even if the estimate for M<j us t depends on the value 
of q, the range spanned by our sample for q = 2.5 and 
q = 3 turns out to be very similaiF^I. namely ~ 2 • 10~ 5 — 
3 • 10 -3 Mq, corresponding to roughly 6 — 1000 M®. By 



15 This discrepancy is probably due to a different choice of the 
literature data used for the derivation of a. In the case of R10 
all the data obtained at sub- mm and mm wavelengths have been 
used. 

16 Note however that for the four disks with the lowest value 
of p (»; 0) Mj~ t cannot be obtained. The reason for this is that 
these very low /3-values cannot be explained with q = 3 (see 
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Table 3. Disk properties. 



Ricci et al.: Dust grain growth in the p-Oph disks 



Object name 


R out (AW07) 


Rout— interval 


Q 


P 




M q=2.5 


M q=3 

iw dust 


M 9=3 - 5 

lw dust 




(AU) 


(AU) 






(M x cm'g- 1 ) 


(M e ) 


(M Q ) 


(M ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 



SR 4 




100-300 


2.5 


0.7 


1.3 


10" 4 


2.1 


10~ 5 


2.4 


10~ 5 


GSS 26 




100-300 


1.9 


0.0 


3.5 


10" 4 


2.3 


10 -3 






EL 20 




100-300 


2.5 


0.8 


3.2 


10" 4 


4.3 


10" 5 


4.5 


10" 5 


DoAr 25 


200 


100-300 


2.3 


0.5 


8.0 


10" 4 


1.8 


10" 4 


2.6 


10" 4 


EL 24 


175 


75-275 


2.2 


0.4 


9.9 


10~ 4 


2.9 


10" 4 


6.3 


10" 4 


EL 27 


275 


175-375 


2.2 


0.5 


1.5 


10" 3 


3.5 


10" 4 


5.6 


10" 4 


SR 21 


600 


500-700 


2.9 


1.1 


5.3 


10" 4 


4.9 


10" 5 


4.5 


10" 5 


IRS 41 




100-300 


2.1 


0.3 


1.3 


10" 4 


6.6 


10" 5 


6.8 


1Q -4 


YLW 16c 




100-300 


2.4 


0.0 


2.3 


10" 4 


4.3 


10" 5 


5.6 


10~ 5 


IRS 49 




100-300 


1.8 


0.0 


3.9 


10" 5 


2.2 


10" 3 






DoAr 33 




100-300 


2.2 


0.4 


1.2 


ur 4 


3.4 


10~ 5 


9.9 


10" 5 


WSB 52 




100-300 


1.8 


0.0 


1.4 


10" 5 


2.6 


10" 3 






WSB 60 


350 


250-450 


1.9 


0.3 


5.6 


10" 4 


2.9 


10" 4 


3.0 


10 -3 


DoAr 44 




100-300 


2.2 


0.4 


3.0 


10" 4 


8.8 


10" 5 


1.9 


10" 4 


RNO 90 




100-300 


2.3 


0.4 


1.1 


10" 4 


3.1 


10" 5 


7.1 


10" 5 


Wa Oph 6 


275 


175-375 


2.4 


0.7 


4.9 


KT 4 


8.0 


10" 5 


9.8 


10" 5 


AS 209 


200 


100-300 


2.4 


0.7 


7.9 


10" 4 


1.2 


10" 4 


1.4 


10" 4 



1.9 ■ 10" 4 



5.5 ■ 10" 5 



1) Underlined objects are those which have not been mapped to date through high-angular resolution imaging. The objects which 
have been mapped but do not have an estimate for the outer disk radius as reported in Column (2) have been spatially resolve by 
Andrews et al. I|2009p . Contrary to Andrews & Williams (12007a| , they modeled the disk surface brightness by using a self-similar 
profile instead of a truncated power-law. For this reason no estimate for 7? ou t could be extracted for these sources (see footnote in 
Sect. 4.1). 2) Best-fit estimate of the disk outer radius by fitting the observed visibilities at sub-millimeter wavelengths using a 
truncated power-law for the surface density profile. AW07: Andrews & Williams (2007a:). 3) Interval of the disk outer radius 
adopted for our analysis. 4) Best-fit estimate of the spectral index of the SED a between 1 and 3 mm derived by considering for 
the outer disk radius the central value of the interval reported in Column (3), and for the power-law index of the surface density 
profile p = 1. 5) Best-fit estimate of the spectral index of the dust opacity ft between 1 and 3 mm. 6) Product between the dust 
mass and the dust opacity at 1 mm from the best-fit two-layer disk model. 7) Dust mass obtained with a power-law index for the 
grain size distribution q = 2.5. 8) Like Column (7) but with q = 3. The sources without an estimate of the dust mass have a 
/3-value (reported in Column (6)) which cannot be reproduced with q = 3. 9) Like Column (8) but with q — 3.5. 



defining a planetesimal as a rocky body with a radius of 
10 km and a density similar to the one adopted for our 
dust grain model (i.e. p ~ 1 g/cm 3 ), the maximum number 
of planetesimals which can be potentially formed out of 
this reservoir of small grains (see discussion below) is ~ 
10 10 -10 12 . 

These numbers have to be taken with great caution 
mainly because of the large uncertainty for the inferred dust 
mass of a factor as large as 10 (see e.g. the discussion in 
Natta et al. |2"0"0"4")) . 

Here it is important to remember that observations at 
(sub-)mm wavelengths are completely insensitive to peb- 
bles/stones much larger than ~ 1 — 10 cm, since the dust 
opacity decreases as the a max of the dust population in- 
creases at sizes larger than the wavelengths of the observa- 
tions. For this reason the dust masses presented here have 



discussion in Sect. 5.1). For the same reason only two disks have 
an estimate for the dust mass with q = 3.5. 



to be interpreted as lower limits for the real total mass in 
solids, since in principle large pebbles/stones or even larger 
rocky bodies like planetesimals may be already present in 
the disk. 

Finally we have investigated relations between dust 
properties in disks (i.e. 0, Md us t) and the stellar ones (listed 
in Table 2) but we did not find any significant correlation, 
similarly to the case of Taurus (R10). No significant cor- 
relation was found neither between grain growth and dust 
mass. 

6. Summary 

We have presented new observations at ~ 3 mm obtained 
with the ATCA array and the new CABB digital filter bank 
for 27 protoplanetary disks in the p— Oph star forming re- 
gion. Among these we selected the 17 isolated class II YSOs 
with well characterized stellar properties (see selection cri- 
teria in Sect. 3.1). Our sample comprises all the p— Oph 
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isolated class II YSOs with an observed flux at 1.3 mm 
larger than ~ 75 mJy, and ~ 50% of the isolated PMS 
stars with mass larger than ~ 0.5 M Q . We have analyzed 
the (sub-)millimeter SED of our disk sample and here is the 
summary of our main findings: 

1. The spectral index f3 of the millimeter dust opacity 
turns out to be lower than the typical value found for 
the ISM for all the 17 disks detected at 3 mm. The 
mean value is < f3 >« 0.5. For the ten disks which 
have been observed and spatially resolved through past 
high-angular resolution continuum imaging at sub-mm 
wavelengths this represents evidence for the presence 
of dust grains as large as at least ~ 1 mm in the disk 
outer regions. For the seven disks which have not been 
mapped yet the observations could in principle be con- 
sistent also with very compact (i? ut < 20 — 30 AU) 
disks, significantly different from those mapped so far. 

2. From a comparison between the results found for our 
sample in p— Oph and an homogeneously selected sam- 
ple of 21 isolated class II disks in Taurus-Auriga (R10), 
there is no statistical evidence of any difference between 
the distribution of /3-values found in the two star form- 
ing regions. This may suggest that environmental ef- 
fects do not play an important role in the first phases 
of planet formation. 

3. There is no evidence for any evolution of the dust spec- 
tral index: dust grains appear to be present in the outer 
regions of protoplanetary disks throughout all the class 
II evolutionary stage of YSOs, confirming what previ- 
ously found in Taurus. In order to explain the retention 
of large dust grains in the outer disk some mechanisms 
which may halt the inward drift of solid particles, e.g. 
local pressure maxima due to turbulent vortices or spi- 
ral density waves, have to be invoked. Since evidence for 
grain growth to millimeter sizes appears to be present 
in some class YSOs the formation of the ~ mm-sized 
grains seen in class II disks may already occur in the 
densest inner regions of molecular cloud cores (Ormel 
et al. 12009)) . 

4. The mm slope of the SED does not correlate with indica- 
tors of dust processing from the silicate feature observed 
at about 10 /im, which are sensitive to grain growth to 
micron sizes in the surface layers of the inner disk; fur- 
ther observations at (sub-) mm wavelengths are needed 
to extend the investigation to a larger sample. 

5. The spanned range in dust mass contained in grains 
with sizes <1 cm as derived with the dust model de- 
scribed in Sect. 4.2 is about 2 • 10~ 5 - 3 • 10~ 3 M G or 
roughly 6 — 1000 M®. This reservoir of small grains is 
capable of forming about 10 10 — 10 12 10 km-sized plan- 
etesimals with a mean density of ~ 1 g/cm 3 . 
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